Symmetry fractional quantization in two dimensions 
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We introduce a solvable spin-rotational and time-reversal invariant spin-1 model in two dimen- 
sions. Depending on parameters, the ground state is an equal-weight superposition of all valence 
loops called "resonating valence loop" (RVL) or an equal-weight superposition of all valence bonds 
known as "resonating valence bond" (RVB). In RVL, ends of open loops are deconfined spinous of 
spin-1/2 that cannot be obtained by simple combinations of spin-1 - a phenomenon of fractionaliza- 
tion; while in RVB, all quasiparticles carry an integer spin, hence no fractionalization. RVL and RVB 
are spin liquids with identical topological order but different spin-rotational and time-reversal sym- 
metry quantum numbers of quasiparticles. We propose that quantized symmetry quantum number 
gives a systematic way to (partially) classify phases with identical topological order in dimensions 
greater than one. 



Quantum number fractionalization [ij is among the 
most striking quantum phenomena in condensed mat- 
ter physics: quasiparticles can have fractional quantum 
number that cannot be obtained by simple combina- 
tion of fundamental constituents of the system. In most 
known examples of dimensions greater than one, frac- 
tional quantum numbers arc quantized to certain discrete 
values and the origin of fractional quantization comes 
from topological order Q of the underlying system. For 
instance, topological order is solely responsible for the 
quantization of fractional charge of quasiparticles in frac- 
tional quantum Hall (FQH) states, whereas the C/(l) 
charge symmetry allows, a priori, continuously tunable 
fractional chargeQ- We call such phenomenon ^topolog- 
ical fractional quantization^^ . 

There are also cases where the symmetry property 
by itself dictates quantization of fractional quantum 
numbers. This occurs when the symmetry group has 
only discrete irreducible (projective) representations. We 
call such phenomenon ^^symmetry fractional quantiza- 
tion" . For both cases of topological and symmetry frac- 
tional quantization, it is possible to realize two dis- 
tinct states exhibiting identical topological order but sup- 
porting quasiparticles with different quantized fractional 
quantum numbers. 

To illustrate symmetry fractional quantization, we in- 
troduce an exactly solvable 50(3) spin-rotational and 
time-reversal invariant spin-1 model on the (decorated) 
honeycomb lattices. Depending on parameters in the 
model, the ground state is an equal-weight superposi- 
tion of valence loops called "resonating valence loop" 
(RVL) or an equal- weight superposition of valence bonds 
known as "resonating valence bond" (RVB)[^, Q. The 
RVL and RVB are Z2 spin liquids with identical topo- 
logical order. But they are distinct phases of matter be- 
cause their quasiparticles carry different quantum num- 
ber quantized by symmetry. Specifically, in RVL, each 



loop is like an AKLT spin-1 chain and the ground 
state is a loop soup pictorially. Ends of open loops are 
deconfined spinous with spin-1/2. Spin-1/2 objects form 
a projective representation of 50(3) group which cannot 
be obtained from any combination (tensor product) of in- 
teger spin representations 0. Moreover, spin-1/2 trans- 
forms differently from spin-1 under time reversal opera- 
tion T: 7"^ = 1 for spin-1, whereas = — 1 for spin-1/2; 
the latter leads to Kramer's degeneracy. Thus the RVL 
exhibits fractionalization of quantum numbers of spin- 
rotational and time-reversal symmetry. To the best of 
our knowledge, this is the first example of quantum num- 
ber fractionalization of time-reversal symmetry realized 
in a solvable modelQ. In contrast, in RVB, all quasipar- 
ticles carry integer spins, hence no fractionalization. 

One motivation to study symmetry fractional quan- 
tization is to classify topological phases of matter with 
symmetry. For free fermion systems, the topological clas- 
sification with discrete symmetries such as time-reversal 
symmetry and particle-hole symmetry has been stud- 
ied systematically [l3l since the recent discovery of topo- 
logical insulators and superconductors [lll - [l3| . More re- 
cently, the classification has been gene ralized to insula- 
tors with certain spatial symmetries [lj-ll6[. For interact- 
ing systems, the classification has been intensively stud- 
ied in lD[iM3; however, classification in higher dimen- 
sions is far from complete (see, nonetheless, Refs. (20l - [23| ). 
We hereby propose that symmetry quantum numbers 
of quasiparticles can (partially) distinguish topologically 
different phases with the same symmetry. The reason is 
that symmetry quantum numbers are discrete quantities 
and cannot change smoothly as long as the underlining 
symmetries are respected. Therefore two gapped phases 
in which quasiparticles transform under different irre- 
ducible representations of a symmetry group cannot be 
adiabatically connected under symmetry-preserving de- 
formations: they are necessarily different phases. Topo- 
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logical order, combined with (fractional) symmetry quan- 
tum numbers, give a more refined classification of topo- 
logi cal p hases with symmetry. It is worth to mention that 
Ref . 22, l23| has classified hundreds of distinct Z2 spin liq- 
uids of spin- 1/2 systems protected by space group sym- 
metry, based on the projective symmetry group (PSG) 
approach to study the variational ground state wave func- 
tion obtained from slave-particle construction. The con- 
nection between the PSG property of ground states and 
the physical quantum number of quasiparticles is not per- 
fectly understood. Moreover, the distinction between 
spin liquids protected by spatial symmetry is not robust 
against disorder, whereas our study dealing with time- 
reversal symmetry is. 

Resonating valence loops (RVL): For simplicity, 
we shall consider a spin-1 system on the honeycomb lat- 
tice first. On each site i, Sf = \b\^a'^^hiv, where cr" are 
Pauli matrices and bi^ are Schwinger bosons with the 



urations: for any loop-covering c, a corresponding dimcr 
configuration is obtained by occupying a virtual dimer 
on any link not covered by the loop; consequently, the 
number of loop configurations is identical to the number 
of dimer ones that increases exponentially with the size 
of the lattice. Since B]- = — Sj^, in Eq. ([T]), we take the 
convention £ A{B) sublattices for all links (ij) such 
that (c| c) < where c is obtained from |c) by flipping 
singlet bonds which are flippable on some hexagon pla- 
quette. The negative overlap 24j between |c) and \c) is 
important for reasons which will become clear later. 

Before dicussing properties of the RVL state, we shall 
first introduce a microscopic model for which the RVL 
state in Eq. ([1]) is the ground state. Since the RVL state 
consists of loop configurations, it is desired to associate 
an energy cost for non-loop states. To achive this, we 
consider the following Klein term[25| 



constraint ^Lfe, 



It^^t + = 2. In the spirit of AKLT, Hk = jJ2 P3W{i)] + PiWii)] + P3[A/'a(j)] 



each of two bosons on each site can form a singlet bond 
with another boson on its nearest neighbor site - each 
spin-1 can participate two singlet bonds. A spin singlet 
bond on the link (ij) is created by Bjj = ^'^'^bl^b^^^ (e^" is 
the Levi-Civita symbol). For reasons which will be clear 
later, we are interested in the following states: 



,(2) 



L LI ± 

where J > 0, A/ (i) labels a cluster of four sites consisting 
of i and three neighbors of i, Ma{i) a cluster of three 
sites consisting of i and two neighbors of i, as shown in 
Fig. [Ija). Here Ps\N'{i)] is a projection operator onto 
the sector of total spin-s of the cluster namely 
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where c is a loop-covering configuration which consists 
of non-intercepting loops and which touches every site 
by one and only one loop and ric counts the number of 
singlet bonds of the loop configuration c on those ver- 
tical links with upper site in sublattice A. |^)r;YL 
an equal weight superposition of all loop-covering config- 
urations, which we propose to call "resonating valence 
loop" (RVL) state. A typical loop-covering configura- 
tion is shown in Fig. [1] For the honeycomb (or any 
other trivalent) lattice, there is a one-to-one correspon- 
dence between loop-covering and dimer-covering config- 
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^ieA/'(i) '^j- PsWa{i)] and other projec- 




FIG. 1: (a) The schematic representation of the honeycomb 
lattice and a typical loop-covering configuration. The black 
(white) sites are in A (B) sublattices. Thick bond represent a 
spin singlet created by Bj^ (see text) and they forms a loop- 
covering on the honeycomb lattice. A hexagon plaquette is 
flippable if it has three thick and three thin bonds. For in- 
stance, the loop configuration is flippable on the hexagon pla- 
quette marked by a. (b) A typecal spin-1 dimer-covering con- 
flguration in which each dimer consists of two singlet bonds. 



where S^^i^ = 

tion operators which will appear later are defined simi- 
larly. Since Hk is a sum of projection operators, it is 
positive semi-definite. It is clear that Hk |c) = 0; |c) 
is a ground state of Hk- As usual, one question about 
the Klein term is whether it is perfect, namely whether 
there are non-loop configurations whose energy is zero. 
Although lacking a rigorous proof 26| , we believe that 
the set of |c) forms the complete ground state manifold 
of Hk since there are no obvious non-loop states with 
zero energy. 

The dual description in terms of dimers tells us how 
loop configurations resonate: if c has three separated sin- 
glet bonds on a hexagon plaquette, c is said to be flip- 
pable on that plaquette, as shown in Fig. [Ija) . It is then 
clear that the maximum possible value of the total spin 
Sa = X^iea ^ hexagon a is 3 when a is flippable. For 
all other configurations non-fiippable on a, the maximum 
possible spin on a is at most 2. This fact makes flippable 
configuration unique and enables us to construct the fol- 
lowing Hamiltonian for which we shall prove that |^')j^vl 
is its ground state. 

^^RVL = Bk + Hq, 

= HK + J'Y.P3{a)[P3{Ea)+P3{Oa)], (4) 

where Oa{Ea) represents the odd (even) sites, namely 
the three A{B) sublattice sites, of plaquette a. Note that 
Pz{a)^ PsiEa) and P3{Oa) commute with each other. 
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It is clear that J^rvl is positive semi-definite since it 
is also a sum of projection operators. Consequently, we 
can conclude that |4')rvl ^ ground state of TJrvl if we 
prove Hrvl |*)rvl = or equivalently TJg |*)pyL =0 
by knowing that Hk |^)j^vl ~ ^- '^'^ prove this, we 
choose a plaquette a and rewrite |'I')pyL = iV'a) + 
J2c I'^o)' "where Ca label loop configurations which are 
flippable on a, and \ipa) contains all other configura- 
tions non-flippable on a. As discussed earlier, Sa < 
2 for a configuration non-flippable on a; this implies 
that P3{a)[P3{Ea) + PsiOa)] IV'a) = 0. For each |ca), 
^3(1) |ca} 7^ 0; however, there is another configuration 
\ca) which is related to \ca) by fiipping the singlet bonds 
on a. Now we shall show that PaiEa) (|ca) -I- |cq)) — 
and the same for PsiOa)- The projection P3{Ea) can be 
written in spin- 3 coherent states 



PsiEa) = ^j<fn I^B^ = 3, n) {Se^ = 3, n| 



(5) 



where \Sec, = 3, n) is the state with Se^ = 3 and 
the spin component S^^ • n = 3 along the unit vec- 
tor n. Define PfiEa) °= \Se, = 3, n) {Se^ = 3, n| as 
the projector to the n direction. Due to spin rota- 
tion invariance of the loop states \ca) and |ca), one 
can prove P^{Ea) (|cq) -I- |ca)) = for aU fi if it's 
true for a particular n. Thus we take n = z and 
study the action of P^{Ea)- From (|ca) + |cq)) = 



B 



Y\.{ij)£Ca.(ij)<^a 

clear that the contribution from 
S^{Ea) = 5| + 5| 
bution from 



to 



Sq is at most 3/2 while the contri- 
16^2^54 is at most 1/2; 
consequently, S^{Ea) < 2 for {\ca) + \ca)), which implies 
that P^{Ea){\ca) + \ca)) ~ 0. Consequently, 
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^fRVL|*)RVL = (^A-+i/Q)|vI') 



RVL 



0. 



(6) 



Thus, we have proved that |^')p{yL is the exact ground 
state of TJrvl. Since P3{a)[P3{Ea) + P^iOaMca) - 
\ca)) 7^ 0, we believe that the ground state of TJrvl in 
each topological sector is unique. Because the honey- 
comb lattice is bipartite, by its analogy to the quantum 



dimer model [27[ the RVL state on the honeycomb lat- 
tice is expected to describe a critical point rather than 
a stable phase - small perturbations would drive it into 
a different state [1^. The exact nature of states at and 
around this critical point remains to be studied in future. 

Decorated honeycomb lattice: To obtain a stable 
spin liquid phase described by the RVL state, we consider 
a spin-1 system on the decorated honeycomb lattice (see 
Fig. 1 in Rcf. [1^), also known as the star or the 3-12 
lattice. The RVL state here is formally the same as in 
Eq. 11]). We consider the following Hamiltonian that is a 
direct extention of the one studied on honeycomb lattice: 

-ffRVL = Hk + J'J2 PeiamiEa) + PeiOa)], (7) 



where Hk is the Klein term which is identical to the 
one in Eq. ([2]). Here a labels a dodecagon instead of 
hexagon plaquette. It is straightfoward to prove that 
^)rvl decorated honeycomb lattice is the ground 

state of i?RVL in Eq. ([7]). Since the decorated honeycomb 
lattice is non-bipartite, by analogy to quantum dimer 
modeljl^l, the short-range RVL state is gapped and then 
stable against small perturbations. We expect that there 
is a finite region in parameter space around the Hamilto- 
nian Eq. (I7|) such that the ground state is a loop-liquid 
state with a finite gap to all excited states. From now on, 
we will focus on the decorated honeycomb lattice unless 
stated othervise. 

Deconfined spinons: In the ground state, all loops 
are closed in each |c). What happens when some loops 
are open? Since each loop is like a spin-1 AKLT chain, 
from the physics of ID AKLT spin-1 chain, it is clear that 
spin- 1/2 excitations would appear at the ends of open 
loops. These spin-1/2 excitations are so-called spinons. 
(In the dual-dimer picture, a spinon excitation occurs at 
the site where two different dimers overlap.) The wave 
function with two static spinons can be written down 
exphcitlyjsij: 



I*, 



biAT. n B, 

Cij {kl)£Cij 



0), 



(8) 



|0), it is where Cy label configurations with only one open loop 
whose ends are i and j. It is clear that l^'i^jv) has two 
spin-1/2 excitations localized at i and j, respectively, and 
each of them is double degenerate. The energy associ- 
ated with the two spinon excitations is finite since the 
constraint imposed by the projection operators is only 
violated around i and j locally evenwhen they are far 
separated. Specifically, the energy cost is at most 4J 
since only P3[Af{i)] and one of PsiJ^faii)] cost finite en- 
ergy around i and similarly around j. Consequently, de- 
confined spinons are supported in the RVL liquid phase 
on the decorated honeycomb lattice. 

Each spinon excitation carries spin-1/2 quantum num- 
ber of the 50(3) spin-rotational symmetry and spin- 
1/2 cannot be obtained from simple combinaitons of 
fundamental constitutents that carry spin-1. Moreover, 
the time-reversal transformation T of spinons satisfies 
= —1. Since the fundamental constituents are spin-1 
which have = 1, spinon excitations exhibit time re- 
versal symmetry quantum number fractionalization. As 
an important consequence, even when the model is per- 
turbed by spin SO(3) non-invariant terms, the double 
degeneracy of "spinon" excitations would remain due to 
the Kramers degeneracy, as long as time-reversal symme- 
try is preserved. 

Resonate valence bonds (RVB): To compare the 
RVL state with other liquid states and obtain better un- 
derstanding on the symmetry fractional quantization, we 
study a spin-1 RVB state which is a direct generalization 
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of the spin-1/2 RVB state given by 

i*)rvb = Ei^)' ic) = (-ir^ n K)'io)' (9) 

c {ij)ec 
where c labels a dimer-covering and Uc is chosen similarly 
such that (c| c) < 0. The spin on each site forms a sin- 
glet with one of its neighbors in this spin-1 RVB state. 
With some straightforward calculation, it is clear that 
this spin-1 RVB state is a ground state of the following 
Hamiltonian with short-range interactions 



H 



RVB 



.s=3 



a = l 



ieOa 



,(10) 



where ii' label the link not in the hexagon plaquette a. 
Again, by analogy with quantum dimer model [s^, we 
believe that this Hamiltonian is gapped and the spin-1 
RVB state is the unique ground state in each topological 
sector. The RVB and RVL states are both Z2 spin liq- 
uids with identical topological order. However, there is a 
qualitative distinction between the two states: spin-1/2 
(spinon) excitations in the spin-1 RVB phase is linearly 
confined while they are deconfined in the RVL phase. 

Phase transition(s) between RVL and RVB: 
When the time-reversal, 5*0(3) spin-rotation, or both 
symmetries are respected, spinous exhibit fractional 
quantization of time reversal or SO(S) quantum num- 
bers; consequently it is impossible to adiabatically con- 
nect the RVL state with the RVB without going through 
a phase transition even though they have identical topo- 
logical order. To study the possible phase transition(s) 
between them, we introduce the Hamiltonian H{X) = 
Ai?RVB + (1 — A)i?R,vL- By tuning A from to 1, we ex- 
pect at least one phase transition between the RVL phase 
at A = and the RVB phase at A = 1. Nonetheless, when 
both time reversal and 5*0(3) are broken, it is expected 
that these two states can be adiabatically connected - 
they belong to the same phase. 

Concluding remarks: In conclusion, RVL and RVB 
states, both of which are Z2 quantum spin liquids, are two 
distinct phases protected by the given symmetry (50(3) 
spin- rotation or time-reversal). We have explicitly re- 
vealed the distinction between RVL and RVB states by 
studying the symmetry quantum numbers of their quasi- 
particlcs. This points to a systematic approach to (par- 
tially) classify symmetric topological phases in dimen- 
sions greater than one provided the generic interplay be- 
tween symmetry quantum numbers and topological order 
(such as fusion rules of quasiparticles) is understood. 
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